Plasmons in graphene: the nearest and next-nearest neighbor hopping 
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The spectrum of plasmon excitations in graphene is studied in the framework of the tight-binding 
model, including nearest-neighbor and next-nearest-neighbor hopping. The dielectric function is nu- 
merically calculated as a function of wave vector and frequency, and two alternative approximations 
for calculating the plasmon dispersion law - the method based on the loss function and that based 
on the nullification of the real part of dielectric function - were used, compared and shown to agree 
with each other. It was also shown that next-nearest neighbor hoping may infiuence substantially 
the plasmon dispersion curve and/or the region of wave vectors where plasmons can be observed. 

PACS numbers: 73.20.Mf 



INTRODUCTION 

Graphene, a single layer of carbon atoms arranged as a 
honeycomb lattice, is a semimetal with remarkable phys- 
ical properties [l|, Q . This is due to the band structure 
of the material which consists of two bands which touch 
each other at two nodes. The electronic spectrum around 
these two nodes is linear and can be approximated by 
Dirac cones. Many unusual physical properties are con- 
trolled by the fact that the Fermi energy is at the nodes, 
where the density of states vanishes. However, graphene 
can be gated, so the Fermi energy can be freely tuned. 

In the following we will discuss collective excitation 
of the 2D electron gas in graphene, caused by an exter- 
nal electromagnetic field. A standard problem in semi- 
conductor physics, it was initially studied in the case of 
graphene only in the Dirac approximation around the 
nodes Later the results of calculations based on 

the full band structure, taken into account within the 
tight-binding model have appeared [1, 0| ■ We extend the 
nearest-neighbor approximation used in these papers by 
taking into account in addition the next-nearest-neighbor 
hopping. 

We consider an electron gas which is subject to an elec- 
tromagnetic potential Vi{q,u!). The response of the elec- 
tron gas is a screening potential Ks(q, uj) which is caused 
by the rearrangement of the electrons due to external 
potential. Therefore, the total potential, acting on the 
electrons, is 

V{q,Lo)^V,{q,u:) + Vsiq,u:). (1) 

Vs can be evaluated self-consistently and is expressed via 
the dielectric function e(q, w). Then the total potential 
reads 

nq,-) = ^^.(q,-). (2) 



ELEMENTARY ELECTRONIC PROPERTIES OF 
GRAPHENE 

Graphene is made out of carbon atoms arranged in 
hexagonal structure. The structure can be seen as a tri- 
angular lattice with a basis of two atoms per unit cell 
(atom A and atom B). The tight-binding Hamiltonian 




FIG. 1: A honeycomb lattice and its Brillouin zone. Left: 
lattice structure of graphene, made out of two 
interpenetrating triangular lattices (ai and a2 are the lattice 
unit vectors, and Si, i=l,2,3 are the nearest-neighbor 
vectors). Right: corresponding Brillouin zone. The Dirac 
cones are located at the K and K' points. 

for electrons in graphene with both nearest- and next- 
nearest-neighbor hopping has the form[lj (we use units 
such that h — 1) 

<i,j>,cr 

-t' E + + i/.c.) (3) 

«i,j>>,cr 

where a„^i{al^ annihilates (creates) an electron with 
spin a{a =t)i) on site Rj on sublattice A (an equiva- 
lent definition is used for sublattice B), t « 2.8 eV is 
the nearest-neighbor hopping energy (hopping between 
different sublattices), and t' is the next nearest-neighbor 
hopping integral (hopping in the same sublattice). The 
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FIG. 2: Energy dispersion of graphene with t' — 0.2t. 
The figure shows a broken particle-hole symmetry. The 
Dirac nodes are shifted up by 3t' . The energy is 
measured in the units of the electronic bandwith 
A = 3t. 



value of t' is not well known but ab initio calculations 
find 0.02i ^ t' ^ 0.2t depending on the tightbinding 
parametrization 

The energy bands derived from this Hamiltonian have 
the form 




i?±(k) =±iv/3 + /(k)-t7(k), 
/(k) — 2 cos{V^kya) + 4 cos(-^fcya) cos(-fcyo) (4) 

where the plus sign applies to the upper (tt) and the 
minus sign the lower (tt*) band. 

The expansion of the spectrum around the Dirac point 
K including t' up to second order in q/K is given by 

E±{k) 3t' ± x;^|q| - ± ^s^ni39^)^ jq^j 

(5) 

where 9q ~ tan^^(2^), and q = K — k. Hence, the 
presence of t' shifts in energy the position of the Dirac 
point and breaks electron-hole symmetry. 



DIRAC POINTS AND SADDLE POINTS 



The Hamiltonian for a general two-sub-lattice lattice 



H 



(6) 



where a is an arbitrary lattice vector, and 5 is a vector 
connecting two sites of different sub-lattices. The disper- 
sion law thus is given by equation 



i?(k)+^t'(a> 



(7) 



i\i.a. 



= 0. 



The conditions for the Dirac points are 
OF 



dE{k) 







OF 



(8) 
(9) 



Actually, in mathematics Dirac points we are dealing 
with are called conical points of the surface, and the con- 
dition dH]) can be found in Ref. Eqs. ^ gives 



(10) 



Substituting this into Eq. ([7]) we get necessary condition 
of the Dirac point 



(11) 



We also see that if k satisfies Eq. (fTT|) . Eq. (|9]) is sat- 
isfied automatically. An interesting fact is that intra- 
sublattice hopping does not influence the position of 
Dirac points, but of course influences the value of en- 
ergy in these points, which in general case is given by 
Eq. (jlO[) : in a particular model we are considering, this 
energy is 3i'. The only symmetry which was essential 
for the Dirac point condition was the sub-lattices 
symmetry (equality of diagonal matrix elements in the 
Hamiltonian ([6]). 

The critical points of the dispersion law surface de- 
scribed by Eq. (|4]) are found from the system 



dm 



= 0, 



9/(k) 

dk,. 



= 0, 



(12) 



which has the solution corresponding to the Dirac points 
and 4 other 



(0,0), 



0, 



2n 
3a 



3a ' 



TT 

3a ' 



The values of / in these points are 
6,-2,-2,-2. 



(13) 



(14) 



The first solution corresponds to maximum of f{k) (mini- 
mum of valence band and maximum of conduction band) , 
the other three solutions are saddle points. The energy 
in these points is 



E = ±t + 2t'. 



(15) 
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DIELECTRIC FUNCTION 

The longitudinal dielectric function was calculated in 
Random Phase Approximation 



e(q,w) = 1 X q,w 

qn 



(16) 



where k is a dielectric constant and x is a polarizability. 
For polarizability we used the Lindhard formula, which 
in our case after some straightforward calculations can 
be reduced to the expression 



x(q, w) = xi (q, ^) + X2 (q, 

with the intraband contribution 
Xi(q,a;) = lmi ^ / i|KkKk+q + l|^ 

, f{<'^E^^) 



and the interband contribution 

X2( 



(q,a;) = lmi ^ / ijK^^Kk+q - 1| 



(17) 



(18) 



(19) 



X- 



(1)^^(2) 



(1) 

k + q 



(2) 



(2) 

k + q 



+ t{ijj - i5) 



where 



Kk = {hi - ih2)/Ek 
t cos 



(20) 



hi 
h2 



{kxd) cos{ V^kyd) + cos{ksd) 
m{kxd) cos{V3kyd) + sin{ksd) 



f(E) = l/(e'^(^"'') + 1) is the Fermi-Dirac distribution 
function, /? — l/ksT, /i is a chemical potential. The 
energies are defined as 



-Ek,i — E, 



(1) 

k,i 



E, 



(2) 
k,i ' 



(21) 



E, 



^k,i 

(2) 
kJ 



4V - + hl= ±tV3T7(k), 



-i7(k) 



('2') (2) 

If we take E^. ' ~ ^k+q = ^ ^^^^^ integral yields the 
same polarizability formula as that found in the nearest- 
neighbor model's polarizability tS, Jj. 



PLASMONS IN GRAPHENE 

Plasmons are known as collective excitations of elec- 
trons, found from the equation 



e(q,^) = 0. 



(22) 



Complex solution of this equation w(q) gives both the 
dispersion and the decay of plasmons. However, numer- 
ical solution of Eq. ([^ for a complicated e(q, w), like 
that given by Eqs. - is challenging. That is 

why approximations schemes are used for analyzing the 
plasmon dispersion. One of them is based upon the loss 
function 



Im 



1 



-/TO[e(q, 



e(q,c.) / {i?e[e(q,c.)]}2 + {/m[e(q, u;)]}^ 



(23) 



calculated as a function of real lu and q. The peaks of the 
loss function are associated with the collective excitations 
of electrons (plasmons), thus giving real function i^(q). 
The other approximation is abridging Eq. ((22)) to the 
equation 



i?e[e(q,c^)] -0, 



(24) 



which is solved for real function w(q). 

In the present paper the polarizability of graphene x 
was evaluated numerically and the corresponding dielec- 
tric function is obtained from Eq. (IT5|) for different values 
of the real frequency oj, the wave vector q and chemical 
potential (Fermy energy) /i (k was taken to be equal to 
4.0). Comparison of our results with those of Ref. Q 
shows that the next nearest-neighbor hopping strongly 
changes the shape of the plasmon curve for the low val- 
ues of chemical potential /i. For high values of chemical 
potentials (higher then 3t', which is equal to the dis- 
placement of the Dirac nodes) the shape of the curves 
does not change much, but the plasmon curve becomes 
shorter. Thus the region of wave vectors, where plasmon 
can be observed, is reduced. 

Plasmons near the saddle points appear on Fig. ^ 
and Fig. ([9]). The results were calculated only for the 
upper energy band where chemical potential is very close 
to the saddle points location. The figures correspond to 
the two values of t' — 0.02, 0.2 and like on other figures 
the direction of qx = 0.0 is selected. The plasmons near 
the saddle points do not have the difference on the shape 
of the curve and look very similar to the other chemical 
potentials of the upper energy band. 

The Fig. (fTO|) compares the results of the two approxi- 
mate methods of calculating the plasmon dispersion law, 
mentioned above. The loss function method has only 
one curve while the real part method has two curves. 
One curve of the real part method closely corresponds 
to that obtained by the loss function method, the other 
does not appear in the case of the loss function method 
calculations. The relevance of each of this curves can be 
analyzed by calculation of the imaginary part of the di- 
electric function in the vicinity of the curve. We see that 
for one curve given by Eq. ([24]) the imaginary part of e 
is very low. And it is the curve which is given by the loss 
function method. In the vicinity of the other curve the 
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FIG. 3: Plasmon dispersion for ^ = 1.44i,t' = 0.02t and 
different values of Qy component {qx = 0) 
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FIG. 5: Plasmon dispersion for /i = —2.82t,t' 
and different values of Qy component {qx = 



= 0.02t 
0) 
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FIG. 4: Plasmon dispersion for fj, = 0.9t, t' = 0.2t and 
different values of qy component {qx = 0) 
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FIG. 6: Plasmon dispersion for ij, = —3.36t, t' = 0.2t 
and different values of qy component {qx = 0) 



imaginary part of e is high, making it irrelevant. Thus 
the loss function is better choice for the calculation of 
plasmons. 
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FIG. 7: Plasmon dispersion for /i — 0.4t, t' = 0.2t and 
different values of qy component (q^ = 0) 
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FIG. 9: Plasmon dispersion for the saddle point 
/i = 0.8t,t' = 0.2t and different values of qy component 
(fe = 0) 
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8: Plasmon dispersion for the saddle point 
0.98t,t' = 0.02t and different values of qy 
component {q^ = 0) 




(a) = -0.94, t' = Omt,qx ■ 
0.1d-i,cZ = 1.42A 




(c) ti = 0.48, t' = 0.02i,g^ 
0.1d-i,d = 1.42A 



(b) il = -1.12, t' = 0.2t,(3a; 

0.1d-i,d = 1.42A 




(d) /X = 0.3, i' = 0.2t, 
0.1d-\(i = 1.421 



FIG. 10: Red curve represents the real part of dielectric 
function and blue curve represents imaginary part. Two 

small windows show areas where the real part of 
dielectric equals to zero. One window(right) shows that 
the imaginary part of dielectric like the real part equal 
to zero. The other window(left) shows that imaginary 
do not equals to zero but real part do. 



